Dynamics of quantum correlations in colored environments 
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We address the dynamics of entanglement and quantum discord for two non interacting qubits 
initially prepared in a maximally entangled state and then subjected to a classical colored noise, i.e. 
coupled with an external environment characterized by a noise spectrum of the form 1//". More 
specifically, we investigate the dynamics for two different configurations of the environment: in the 
first case the noise spectrum is due to the interaction of each qubit with a single bistable fluctuator 
with an undetermined switching rate, whereas in the second case we consider a collection of classical 
fluctuators with fixed switching rates. Since environmental noise is introduced by means of stochastic 
time-dependent terms in the Hamiltonian, we are able to describe the effects of both separate and 
common environments. The full quantum dynamics of the system is obtained by averaging over 
different noise parameters. We show that environments with the same power spectrum, but different 
configurations, give raise to opposite behavior for the quantum correlations. In particular, depending 
on the characteristics of the environmental noise considered, both entanglement and discord display 
either a monotonic decay or the phenomena of sudden death and revivals. Our results show that 
the microscopic structure of environment, besides its noise spectrum, is relevant for the dynamics of 
quantum correlations, and may be a valid starting point for the engineering of colored environments. 

PACS numbers: 03.65.Yz,03.67.-a,05.40.-a 



I. INTRODUCTION 



Entanglement and discord describe remarkable fea- 
tures of quantum systems. Indeed, they are closely re- 
lated to the amount of quantum correlations contained 
in the system, coming either from non-separability or the 
impossibility of local discrimination Beyond their 
specific role in fundamental physics, entanglement and 
discord has been also recognized as resources for quantum 
technology, e.g. for the processing of quantum informa- 
tion and for the effective implementation of quantum en- 
hanced protocols. In particular, in the last decade it has 
been recognized that separable mixed states may 

represent a resource, if they show a nonzero quantum dis- 
cord. In fact, mixed separable states with non-zero dis- 
cord have been exploited to achieve a speed-up for certain 
computational tasks compared to classical states 0] ■ 
This is true also for continuous variable systems, where 
Gaussian quantum discord has been introduced Q, 
measured 10 
tocols [nl 12 



and exploited for quantum enhanced pro- 



An essential ingredient to exploit the quantumness of 
a physical system is the preservation of its coherent time 
evolution. On the other hand, the unavoidable interac- 
tion with its environment usually destroys coherence and 
quantumness [ij, [13], and in turn its use for quantum 
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technology. For these reasons, much attention has been 
devoted to the analysis, characterization and control of 
the dynamics of quantum correlations in different phys- 
ical systems, including quantum optics [H, [l6, 18, 1^, 
nuclear magnetic resonance [2l| , nanophysics (2^ ^J] 
and biology [13, H^. 

For bipartite open quantum systems interacting with 
a quantum environment, entanglement and discord may 
exhibit peculiar features such as sudden death and tran- 
sitions, revivals, and trapping [26- 31]. Such phenomena 
have been linked either to direct [32] or indirect [33] ef- 
fective two-qubit interactions. For non-interacting qubits 
they are due to the non-Markovian nature of the envi- 
ronment ^3^] , which results in the transfer of correlations 
back and forth from the two-qubit system to the various 
parts of the total system. Recently, revivals of quantum 
correlations have been found also for quantum system 
coupled to classical sources jH, and have been con- 
nected to a quantifier of non-Markovianity for the dy- 
namics of a single-qubit. Indeed, it was proven that a 
classical noise can mimic, without loss of generality, a 
quantum environment not affected by the system or in- 
fiuenced in a way that does not result in back-action [s^ . 

Among the class of open quantum systems interact- 
ing with a classical environment, a particular attention 
has been devoted to systems made of two qubits sub- 
jected to a classical source of random telegraph noise 
(RTN) [36l - l4H . namely interacting with a bistable fluc- 
tuator, randomly switching between its two states with 
a given rate 7. Depending upon the ratio between the 
switching rate and the system-environment coupling, the 
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dynamics of a quantum system may exhibit Markovian 
or non-Mar kovian behavior. In fact, the great inter- 
est in RTN is due to the fact that it is able to model 
environmental fluctuations appearing in many nanode- 
vices based on semiconductors, metals and superconduc- 
tors [13, H^. Furthermore, it also represents the basic 
building block to describe noises of the type 1//" , which 
are responsible for decoherence in quantum solid-state 
devices [iil - lisj . This kind of noise spectra stem from 
the collection of random telegraph source with different 
switching rates and are usually referred to as colored spec- 
tra. The color of the noise depends upon the value of the 
parameter a [1^ . For a — 1 the so-called pink 1// noise 
is found, which is obtained from a set of random tele- 
graph fluctuators weighted by the inverse of the switch- 
ing rate. Another interesting case is the 1//^ spectrum, 
also called brown noise from its relation to a Brownian 
motion. 

Environments characterized by a noise spectra usu- 
ally arise when a system is coupled to a large number of 
bistable fluctuators, with a specific distribution of their 
switching rates. Upon considering a collection of fiuctua- 
tors, the colored noise may be implemented by means of 
a linear combination of sources of RTN, each character- 
ized by a specific switching rate chosen from a suitable 
distribution. On the other hand, in this paper we show 
that the same spectrum can be obtained if we consider 
a single fiuctuator with a random switching rate. More- 
over, at variance with the general belief we show that the 
two configurations lead to different physical phenomena, 
i.e. even if the spectra has the same 1//" behavior, the 
dynamical evolution of quantum correlations may be very 
different. In particular, we investigate the time evolution 
of entanglement and discord of two initially entangled 
non-interacting qubits coupled to a classical environment 
described either by means of a single random fiuctuator 
or a collection of RTN sources. For the model of a single 
fiuctuator with random switching rate, the correlations 
decay with a damped oscillating behavior. In the config- 
uration with many bistable fiuctuators, pink noise leads 
to a monotonic decay of entanglement and discord, while 
the presence of brown noise induces phenomena of sud- 
den death and revivals. We ascribe this discrepancies to 
the different number of decoherence channels in the two 
configurations. In other words, the time-evolution of the 
system is determined not only by the spectrum of the 
environment, but also by its configurations, i.e. to its 
microscopic structure. 

For both the configurations, the dynamics of the two 
qubits is ruled by a stochastic Hamiltonian with time de- 
pendent coupling. The average of the time-evolved states 
over the switching parameters describes the dynamics of 
the two qubits. Upon a suitable choice of stochastic time- 
dependent terms in the Hamiltonian, we are able to de- 
scribe the effects of both separate and common environ- 
ments: in the former each qubit is locally coupled to a 
random external signal, while in the latter both qubits 
are subject to the effect of a common classical environ- 



ment. 

The paper is organized as follows: in Section 2, we 
briefly review the definition of negativity as a measure 
of entanglement, and of quantum discord. In Section 3 
we present the physical model for two qubits interacting 
with a classical environment, being represented either by 
a single random bistable fiuctuator or a collection of fiuc- 
tuators. Section 4 reports results for the different con- 
figurations, whereas Section 5 closes the paper with a 
discussion and some concluding remarks. 



II. ENTANGLEMENT AND QUANTUM 
DISCORD 



In this section we briefiy review the concepts of en- 
tanglement and quantum discord between two qubits. 
In particular we evaluate entanglement by means of 
negativity [10], which is given by: 



N ^2 



(1) 



where A~ are the negative eigenvalues of the partial 
transpose of the system density matrix. Note that the 
negativity is bound between zero, for separable states, 
and one, for maximally entangled states. The concept 
of negativity has recently been extended to the case of 
tripartite systems 
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The quantumness of a system is described by discord 
d, [sll , namely the difference between the total and the 
classical correlations in a system. The quantum mutual 
information quantifies the total correlations in a system 
and is defined as: / = S{p^) + S{p^) — S{p) where 
pA{B) jg ^YiQ partial trace of the total bipartite system p 
and S{p) — — Tr/5log2(p) is the von Neumann entropy. 
The classical correlations are evaluated by means of 
the expression [52] C = max{n,} [S{p^) - S[p'^\{Wj})\ 
where {IIj} are projective measurements on subsys- 
tem B and S{p^\{Pij}) = EjPjPf, with pf = 
TrBpjpIIjJ/TrpjpIIj] is the remaining state of A af- 
ter obtaining outcome j on B. Therefore the quantum 
discord is the difference between the mutual information 
and the classical correlation: 



Q^I-C. 



(2) 



Usually, to compute the quantum discord is not an easy 
task, since it involves an optimization problem. But in 
the case of two-qubit X-matrices that can be written as 

p = J (l + J2^=i '^j'^j ® '^j) ' where aj are the three Pauli 
matrices, an analytic expression for Q was derived by Luo 
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Q - Ci - C2 - C3) l0g2(l - Ci - C2 - C3) 

+ (1 - Ci + C2 + C3) l0g2(l -C1+C2+ C3) 
+ (1 + Ci - C2 + C3) l0g2(l +C1-C2+ C3) 
+ (1 + Ci + C2 - C3) log2(l +C1+C2- C3)] 

-^I0g2(l-c)-^l0g2(l+c), (3) 

where c := niax{|ci|, |c2|, |c3|}. 



III. THE PHYSICAL MODEL 

We consider the quantum correlations, in the form 
of both entanglement and discord, between two non- 
interacting qubits subject to noisy environments with 
different configurations. The two qubits are initially pre- 
pared in the Beh state |</)+) = (|00) + |ll))/\/2.' The 
interaction between the qubits and the environment can 
be either local or global, i.e. we examine both the case 
of independent environments acting locally on each qubit 
and of a common environment affecting the two qubits. If 
we set S = 1 and adopt the spin notation, the two-qubit 
Hamiltonian reads: 



H{t)^HA{t)' 



(4) 



where Hj^(^g^ is the Hamiltonian of a single qubit sub- 
ject to a classical time-dependent noise which affects the 
transition amplitude parameter CA{B){i)'- 



HA(B)it) = elA{B} + VCA(B){t)<^x A(B), 



(5) 



with e the qubit energy in absence of noise (energy de- 
generacy is assumed), ^a(b) the identity matrix for sub- 
space A{B), V is the coupling constant between the sys- 
tem and the environment, Gx the Pauli matrix. If the 
time-dependent coefficient c^(B)(i) can randomly flip be- 
tween two values c(t) = ±1 with a fixed rate 7, then 
Eq. (E) d escribes a qubit subject to a random telegraph 
noise^l 113, S El, 111 [Eg. This Hamiltonian, ex- 
tended to describe qudits, has also been used to analyze 
the time evolution of entanglement of a continuous-time 
quantum walk of two indistinguishable particles on a one- 
dimensional lattice [5^ . The Hamiltonian (jlj is stochas- 
tic due to the random nature of the noise parameter c(t) . 
For a specific choice of c(t), the total system evolves ac- 
cording to the evolution operator e~*-^^'* ' with pos- 
itivity ensured by the very structure of the Hamiltonian 
in Eq. ([S]) (s^]. By averaging the global state over dif- 
ferent realizations of the sequences of c(t) , the two-qubit 
mixed state is obtained. 

In order to reproduce the 1//" spectrum, the single 
RTN frequency power density must be integrated over 
the switching rates 7 with a proper distribution: 



(6) 



where S^tn (/, 7) is the random telegraph noise frequency 
spectral density with Lorentzian form S^tn (/, 7) = 
47/(47r^/^ -)- 7^). The integration is performed between 
a minimum and a maximum value of the switching rates, 
respectively 71 and 72. ^(7) is the switching rate distri- 
bution and takes a different form depending on the kind 
of noise: 



7 111(72/71) 



Pa (7) 



(7172)° 



_ T ■sr^ 

72 -7i 



1< a ^ 2 



(7) 



It follows that, in order to simulate a frequency spec- 
trum proportional to 1//", the switching rates must be 
selected from a distribution proportional to I/7". When 
the integration in Eq. ([5]) is performed, the spectrum 
has the requested 1//" behavior in a frequency inter- 
val, so that every frequency belonging to such interval 
satisfies the condition 71 ^ / ^ 72- Eq. ^ can be 
obtained either considering a single bistable fluctuator 
whose switching rate is randomly chosen from a distri- 
bution Pq, (7) or from a collection of sources of RTN each 
with a switching rate taken from the same 7-distribution. 
Even if the spectrum is the same, the physical systems 
described are indeed very different. 



A. 1//" noise from a single fluctuator 

In the case of a single fluctuator, the noise parameter 
c(t) can only flip between two values ±1. The difference 
with the RTN case, is that here the switching rate is not 
known a-priori. This means that the bistable fluctuator is 
described by a statistical mixture whose elements are cho- 
sen from the ensemble {7,^0(7)}- This model describes 
the physical situation in which each of the two qubits 
interacts with a single fluctuator, for separate baths, or 
with the same fluctuator in the case of a common bath. 
The qubits are affected only by one source of noise, and 
therefore only one decoherence channel is present. 

The global system evolves according to the Hamilto- 
nian (U) , with a specific choice of both the parameter c(t) 
and of its switching rate. In particular, for each value of 
the switching rate picked from a distribution ^0(7), we 
create a large number of sequences c(t) . This means that 
we are actually introducing some uncertainty on the rate 
7. This uncertainty gives rise to a 1//" spectrum. For 
each selected switching rate, the evolution corresponds 
to the one of a bistable fluctuator with a characteristic 
RTN phase shown to be: 



<PA(B)(<) ^ \ c{t')dt' 
Jo 



(8) 
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and characterized by a distribution [5J|: 



ten as: 



{ [S{ip + vt) + 5{ip - vt)] + 'l[<d{,p + vt) + Q{>p- vt)] } 
1 (7V1 - i.'P/^iY 



(9) 



where Iv{x) is the modified Bessel function and 9(x) is 
the Heaviside step function. 

For a given 7 the global system is described by the X- 
shaped density matrix obtained by averaging over the 
noise phase the density matrix p{ip, 7, t) corresponding 
to a specific choice of the parameter c(t) [s^ : 



(10) 



Pilit) = J p{if,-f,t)p{ip,t)dip ■ 

^[(l+/3de(ce))|0+>(<^+| 
+ (l-/3de(ce))IV'+)(V'+l 



(11) 



where = (|01) + |10))/\/2 is a Beh state and Pde = 
D2y{t) and Pee — D4^(t) are time-dependent coefficients. 
The function Dmi^{t) represents the average of the RTN 
phase factor i.e. 

(e™^(*)) = / e""^(*)p(^, t) = D^,{t), (12) 



where 



-It 



-It 



cosh (nmut) + sinh (Hmut) 

cos {Km J) + ^ sin {Km J) 



(13) 

for 7 > m;^ and 7 < mi' respectively, where 
Kmi'=\/\'y^ — (mi^)2| with m g {2,4}. The two-qubit 
density matrix is obtained by averaging the density op- 
erator in Eq. ()lip over the switching rates: 



Pde{t) 



Pde{l,i) Pa{lA)Pa{lB)A-tA'^lB (14) 



71 •'71 
72 /•72 



Pee 



(15) 



71 ■'71 



Once the average over 7 is performed, the time-evolved 
density matrix reads: 



1 



Pde(ee)W = 3 [(1 + A<ie(ee) ) I ^ ) I 

+ (l-Arfe(ee))IV^+)(V^' 



(16) 



(■72 

Ade = / D2u Pa (7) (^7 
-•'71 
/•72 

Aee = / DiyPa{'^) 
J-fi 



(17) 



This means that the quantum system is again described 
by an X-form density matrix. 



B. 1//" noise from a collection of fluctuators 

The 1//" noise spectrum can arise from the coupling of 
a system with a large number of fluctuators, each char- 
acterized by a specific switching rate, picked from the 
distribution (7) [13, HH in a range [71,72]. In this case 
the random parameters in Eq.([5]) describes a linear com- 
bination of bistable fluctuators c{t) = J2j=^i where 
Nf is the number of fluctuators and we drop the sub- 
script A{B) to simplify the notation. Each Cj{t) has a 
Lorentzian power spectrum whose sum gives the power 
spectrum of the noise: 



Nf 



In fact, the sum X]j2i 'S'j (/; 7j), with the 7^ belonging 
to the distribution pa{'yj), can be viewed as the Monte- 
Carlo sampling of the the integral Nf J^^ S{f )pa{j) d'j, 
which is, but for a constant, the one in Eq. There- 
fore, in order to obtain a 1//" spectrum, it is necessary 
that a sufficiently large number of fluctuators is consid- 
ered, and that the selected jj are a representative sam- 
ple of the distribution Pa{jj) in the range [71,72]- This 
means that the minimum number of fluctuators we can 
sum- up depends on the range of integration. In fact, a 
big number of fluctuators is required to sample the distri- 
bution of the switching rates over a large range, while few 
fluctuators are sufficient in the case of a narrow range. 
Note that we assume that all the fluctuators have the 
same coupling constant with the environment, that is 
= V for j = 1 . . . N f. 

The global evolution operator U{t) cx e^'^j for 
fixed values of the parameters associated to each fluc- 
tuator, permits to compute the density matrix of the 
global system as a function of a total noise phase Lp(t) — 
'Y^- ifij {t) . Following the approach used in [55| to evalu- 
ate the dynamics of two qubits subject to a single RTN, 
the time-evolved density matrix of the system at time t 
can be expressed as: 



i 



t7|+^47r2/°'/"- (18) 



p{t) 



p{(p,t)pT{ip,t)dip, 



(19) 



where the time-dependent coefficient A^eice) can be writ- 



where pT{<f,t) = Y\jP{>fij,t) is the global noise phase 
distribution. The density matrix in Eq. (IT^ depends 
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on the average of the phase factor ^ which can be 

computed in terms of the RTN coefficient D{t) as: ([T^. 
We have 



(20) 



where the last equahty holds since the RTN phase coeffi- 
cients are independent. By inserting Ea. (l20)) in Eg. dTOl) . 
we evaluate the two-qubit density matrix: 



Pde{ce){'t, {7j}) 



1 r 

2 



(l+r<ie(ce))|0+>(0+| 
+ (l-rrfe(ce))|^ + )(^^ 

The coefficients appearing in the density matrix are 

Tde^YlDjADjBit) with D{t)^D2u{t) 

j 

Fee = ]^i:)j(<) with D{t) = Di^{t) . 



(21) 



(22) 



Note that the fluctuators have fixed switching rates {7^} 
rates {7^}, j = l...Nf. 

For the sake of completeness, here we report also a third 
scenario which describes an environment with a 1//" 
spectrum. This is the case of two qubits subject to a 
collection of fluctuators whose switching rates are known 
with an uncertainty determined by the I/7" distribution 
in the finite range [71 , 72] . The sources of RTN do not 
have a fixed 7, so they are described by an ensemble of 
many fluctuators. The dynamics is evaluated by averag- 
ing the two-qubit density matrix for specific values of the 
switching rates (j2ip over the 7 in a range [71,72]: 



Pde{ce) (*) 



Pde(ce)(i, {7j})Pa({7j})rf{7j}, (23) 



where Pa{{lj}) = YljPailj) and ^{7^} = Ilj /«^7j- The 
density matrix preserves again the X form during time 
evolution. It has the same functional expression of Eq. 
(PT|) . but with coefficients 



r, 



de{ce) 



^'de(ce) — [Ade(ce)] 



Nf 



(24) 



where Kfi^^^ce) are given in Eq. (jl7p . 



range [71,72] of the switching rate. Following the defini- 
tion of negativity in Eq. ([T]), the negativity reads: 



Ndeit) = Adeit) 
N,,{t) = \Aceit)\, 



(25) 
(26) 



where the As are given in Eq. (jl7l) . Since the density 
matrix preserves its X form during time evolution, the 
discord can be evaluated using Eq. Q: 



where 



h{x) 



}de(ce){t) — h{^de{ce)) i (27) 



(1 + X) log2 (1 + X) + (1 - X) l0g2 (1 - X) 



The dynamics of N and D is shown in Fig [T] for the 
cases of pink and brown noise. The range of integration is 
[10^4, 10^]/!/. Negativity and discord have been obtained 
analytically, except for the integral A(je(ce): which must 
be computed numerically. 




FIG. 1: Top panels: Time evolution of negativity (left) 
and discord (right) for two qubits subject to a random 
bistable fluctuator with spectrum 1/f for different (solid line) 
and common (dashed line) environments when [71,72]/;^ ~ 
[10"'*, 10'']. Bottom panels: Time evolution of negativity 
(left) and discord (right) for two qubits subject to a ran- 
dom fluctuator with spectrum for difi^erent (solid line) 
and common (dashed line) environments when [71, 72]/*' = 
[10-^10*]. 



IV. RESULTS 

In this section we present the analytic expression for 
negativity and discord in both the cases of qubit inter- 
acting with a single random fluctuator and a collection 
of Nf bistable fluctuators. 

A. Single fluctuator 

In the case of two qubits interacting with a single ran- 
dom fluctuator, the dynamics depends upon the selected 



Quantum correlations decay with damped oscillations. 
In the case of a qubit subject to different environments 
with 1/f spectrum, the oscillations have a periodicity of 
TT. This periodicity can be explained by analyzing the 
analytic expressions of the quantum correlations. In par- 
ticular, we note that in the integral of Eq. (fT7|) . the 
functions exhibit damped oscillations for 7 < 2z/ with 
periodicity 2Tr/ Kmu, and for 7 > 2;/ monotonically decay. 
Their weighted superposition leads to an interference ef- 
fects that can be summarized as follows: the oscillating 
components result in the formation of alternatively pos- 
itive and negative peaks spaced by 7r/2. On the other 
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hand, the monotonic decay components combine to can- 
cel the negative peaks and preserving the positive ones. 
Finally we are left with an oscillating function with pe- 
riodicity of TT. The same concept applies in the case of 
a common environment, but now the D^i, sum up in an 
oscillating function with periodicity of 7r/2. 

The 1//^ noise spectrum leads to oscillating functions 
of time with periodicity 7r/2 and 7r/4 respectively for dif- 
ferent and common environments and again this peri- 
odicity is related to the fact that with such a distribu- 
tion, the selected values of 7 accumulate near the lower 
value of the frequency range, thus leading to a beat phe- 
nomenon with constructive interference with the above 
mentioned periodicity. If different ranges of integration 
are considered, different time-behavior for the quantum 
correlations can arise, but we will not analyze this effects 
in this paper. 

B. Collection of fluctuators 

In the case of two qubits interacting with a collection 
of fluctuators with fixed switching rates, the dynamics is 
very different, depending on the spectrum of the noise. 
The analytic expression for the negativity and the discord 
is computed starting from the density matrix in Eq. (|2ip . 
The negativity reads: 

Nde{t) = Tdeit) 

iVee(t) = iFceWI, (28) 

where the F's are given in Eq. (|22|) . Also in this case the 
density matrix preserves the X form during time evolu- 
tion, such that the discord is calculated by applying the 
Luo formula ([3]): 

Qde{ce){t) — H^deice)) ■ (29) 

The negativity is the product of many oscillating coeffi- 
cients, with various periodicities. In the case of switching 
rates taken from a I/7 distribution, the product of these 
terms results in a monotonic decay for both entangle- 
ment and discord. In Fig. [5] we report the behavior of 
such quantities in the case of 20 and 100 fluctuators. As 
the number of fluctuators is increased the quantum cor- 
relations decay faster. We consider 20 sources of RTN 
as the minimum number of fluctuator needed to obtain 
both a reliable profile of the frequency spectrum and a 
representative sample of the ^(7) distribution. Though 
it is possible to obtain a pink noise spectrum even with 
a smaller number of fluctuators, this approximation does 
not describe a sample of l/7-distributed switching rates. 
A very different behavior arises when the 7's are selected 
from a 1/7^ distribution. Phenomena of sudden death 
and revivals appear for both entanglement and discord. 
As the number of fluctuators is increased, the heights of 
the peaks decrease. The peaks have a periodicity of 7r/2 
and 7r/4 for different and common environments respec- 
tively. As in the case of the single random fluctuator. 




FIG. 2: Top panels: Time evolution of negativity (left) and 
discord (right) for two qubits subject to a collection of Nf 
bistable fluctuators with spectrum 1/f for different environ- 
ments when [71,72]/;^ = [10^*, 10^]. Bottom panels: Same as 
before but with qubits subject to a common environment. 



this is explained by considering that the selected switch- 
ing rates have small and very close values. The product 
of functions with almost the same periodicity gives a pe- 
riodic behavior with narrow peaks. 

Our results clearly show that the sheer knowledge of the 
spectrum is not sufficient to determine the dynamical 
evolution of correlations. Indeed, it is also the the num- 
ber of decoherence channels that plays a key role. Dif- 
ferent physical models of environments can lead to the 
same spectrum. But, if the two-qubit system interacts 
with only one decoherence channel, then revivals appear 
because the system is affected only by one source of clas- 
sical noise and the information can flow back. If many 
sources of decoherence are present, then the information 
can be completely lost, depending on the channel char- 
acteristics, that is the distribution of the switching rates. 

For the sake of completeness, we report also the ex- 
pressions for negativity and discord in the case of two 
qubits subject to a collection of bistable fluctuators with 
stochastic switching rates. From Eq. (|24l) . we can write 
the negativity and discord as: 

N,,it)^\rUt)\ = \[Adeice)ff\ (31) 
Q<ie(ce)(0=Mr:Ze(ce))> (32) 

The quantum correlations decay exponentially in the case 
of 1// noise and with smooth damped oscillations in the 
case of 1//^ noise. The latter is the most general case, 
in which the collection of fluctuators is described by a 
statistical ensemble. Also in this scenario, brown noise 
leads to non-monotonic decay of negativity and discord, 
since all the selected Dmu{t) terms have similar periodic- 
ity and beat effects arise, with constructive interference. 
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FIG. 3: Top panels: Time evolution of negativity (left) and 
discord (right) for two qubits subject to a collection of TV/ 
bistable fluctuators with spectrum 1/ for independent envi- 
ronments when [71, 72]/!^ = [10~*,10'*]. Bottom panels: Same 
as before, but with qubits subject to a common environment. 



V. CONCLUSIONS 

We have analyzed in details the dynamics of entan- 
glement and discord of a two-qubit system initially pre- 
pared in a maximally entangled state and then subjected 
to classical noise. In particular, we have addressed envi- 
ronments characterized by a 1//" noise spectrum, since 
it is ubiquitous in solid state devices. We analyzed both 
the case of independent environments associated to each 
qubit and of a common environment, and have take into 
account explicitly the structure of the environment. More 
specifically, we have analyzed two distinct configuration 
of the environment, both characterized by a 1//" spec- 
trum, but with a different number of decoherence chan- 
nels. Analytic expressions have been found both for the 
negativity and the discord, which show the same qual- 
itative behavior in all the analyzed scenarios. Indeed, 
discord is a function of negativity only. 

In the first configuration, the two qubits interact with a 
single bistable fluctuator, which has a random switching 
rate leading to an overall 1//" spectrum. An oscillat- 
ing behavior of quantum correlations for both pink and 
brown noise is found. The effect of a common environ- 
ment is to better preserve the correlations and to double 
the number of revivals compared to the case of indepen- 
dent baths. In the presence of pink noise, the heights of 
the peaks of negativity and discord decay faster than in 



the presence of brown noise. In fact, the switching rates 
are more evenly distributed over the range [71,72] than 
in the case of the 1//^ noise, thus leading to massive 
destructive interference. In the second configuration the 
two qubits interacting with a collection of bistable fluc- 
tuators, each one with fixed switching rate chosen from a 
distribution 1 /7" ; also in this case the overall noise spec- 
trum is 1//". In the case of pink noise, entanglement and 
discord show a monotonic decay, while for brown noise, 
sudden death and revivals occur. Quantum correlations 
can be written as a product of oscillating and exponen- 
tial functions. Since the switching rates of the fiuctua- 
tors are selected from a distribution 1//", for the pink 
noise they lead to destructive interference, while brown 
noise enhances constructive interference. As expected, 
the action of independent or common environments has 
different effects on the robustness of correlations in agree- 
ment with previous results obtained with different noise 
models [28]. 

Our results clearly show that the behavior of quantum 
correlations in not infiuenced only by the spectrum of 
the environment, but also by the number of decoherence 
channels, i.e. by the very structure of the environment. 
With a single decoherence channel, revivals of correla- 
tions appear, indicating that a back-fiow of information 
from the environment to the system is possible. When 
the number of decoherence channels is increased, the in- 
formation is quickly lost, and no revivals can occur. 

The characterization and the control of quantum corre- 
lations are fundamental for the development of quantum 
technology. Not only quantum correlations constitute a 
resource to process quantum information, but the deep 
understanding of their nature will provide a better in- 
sight on the nature of quantum states themselves and 
on the transition between quantum and classical descrip- 
tion of physical systems. Our results indicate that the 
microscopic structure of environment, besides its noise 
spectrum, is relevant for the dynamics of quantum cor- 
relations, and may be a valid starting point for the engi- 
neering of colored environments. 
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